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Abstract 

It is known that classical solutions to the one-dimensional quasilinear Smoluchowski-Poisson 
system with nonlinear diffusion a{u) = (1 + n)~^ may blow up in finite time if p > 1 and exist 
globally if p < 1. The case p = 1 thus appears to be critical but it turns out that all solutions 
are global also in that case. Two classes of diffusion coefficients are actually identified in this 
paper, one for which all solutions to the corresponding quasilinear Smoluchowski-Poisson system 
QQ ■ are global and the other one leading to finite time blow-up for sufficiently concentrated initial 

' data. The cornerstone of the proof are an alternative formulation of the Smoluchowski-Poisson 

m ■ 

' system which relies on a novel change of variables and a virial identity. 

, 

1 Introduction. 

: 

I In one space dimension, the quasilinear Smoluchowski-Poisson (SP) system reads 

; dtU = d-, {a{u) d^u - u d^v) in (0,T)x(0, 1), (1) 

= d^v + u-M in (0, T) x (0, 1), (2) 

d,A{u){t,0) = d,Aiu)it,l) = d,vit,0)=dMtA) = for te(0,T) (3) 

u{0) = no > in (0, 1), {v{t)) = for t G (0,T), (4) 



*Institutc of Applied Mathematics, Warsaw University, Banacha 2, 02-097 Warszawa, Poland. E-mail: 
T.Cieslak@impan.gov.pl 

tJnstitut de Mathcmatiqucs de Toulouse, CNRS UMR 5219, Universitc dc Toulouse, F-31062 Toulouse ccdex 9, 
France. E-mail: Philippe.Laurencot@math.univ-toulouse.fr 



1 



where (m, f) : (0,T) x (0, 1) — M are the unknown functions, a G C^((0, oo)) is a positive function, 
A is an indefinite integral of a, and the constant M := (uq), where (/) := Jj f{x) dx/\I\ denotes the 
mean value of a given function / G L^{I) defined on an interval / of M of length |/|. Such systems 
appear in the modelling of self-gravitating systems in astrophysics [5] and chemotactic processes in 
biology [12] . One particular feature of this system is that there is a competition between the diffusion 
dx{a{u)dxu) and the drift term dxiud^v) that may result in the occurrence of blow-up in finite time, a 
phenomenon related to the gravitational/chemotactic collapse in astrophysics and biology. Roughly 
speaking, it has been shown in any space dimension d > 1 that, if a(r) < C (1 + r)~^ for r > and 
some p > {2 — d)/d, there are solutions to the quasilinear SP system blowing up in finite time while 
all solutions are global if a(r) > C (1 + r)~^ for r > and some p < {2 — d)/d ^ (see also [6l [18] for 
a related system in Mf^). Thus, the exponent p = {2 — d)/d appears to be a critical exponent as it 
corresponds to the borderline growth for the nonlinear diffusion coefficient that separates two different 
behaviours. Furthermore, ii d > 2, it has been shown in (TJ [131 112] that, if a(r) = (1 + r)'-'^"^-*/'^, 
there is a threshold value of M above which finite time blow-up takes place for some initial data 
while global existence is guaranteed if M is sufficiently small (see also [2], [6l [TTl [18] and the references 
therein for the quasilinear SP system in Mf^). The one-dimensional case d = 1 (corresponding to 
a(r) = 1/(1 + r)) is not covered by the above mentioned results and the main purpose of this paper 
is to elucidate what happens in that particular case. 

As we shall see below, it turns out that the situation for a{r) = 1/(1 + r) in one space dimension 
is strikingly different from that encountered in higher space dimensions and that all solutions are 
global. We actually identify a class of diffusion coefficients encompassing the one found in [9] for 
which all solutions to ([I])-(|ll) are global. Our approach relies on a new formulation of the SP system 
([T])-(|ll) which transforms it into a single quasilinear parabolic equation. While it is well-known that 
the indefinite integral (with respect to space) of u vanishing at zero solves a quasilinear parabolic 
equation [9l [H] , we go one step beyond and derive the equation satisfied by the first derivative of 
the inverse of the indefinite integral of u. The latter is the solution of a nonlinear reaction-diffusion 
equation with homogeneous Neumann boundary conditions, the reaction term being affine and the 
diffusion term nonlinear. This transformation is described in details in Section [2l Let us emphasize 
at this point that the blow-up in finite time of a solution to ([I])-(|4]) amounts to the so-called "touch- 
down" of the solution to the transformed equation: in other words, a solution to ([I])-(|l]) blows up 
in finite time if and only if the corresponding solution to the transformed equation reaches zero in 
finite time. Section [3] is then devoted to the local-in-time existence and uniqueness of solutions. At 
this point we take advantage of the transformation introduced in Section [2] to provide an simple 
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proof of the uniqueness which is somehow related to the recent use of Wasserstein distances to study- 
nonlinear diffusion equations [3llll[T0]. In Section |U we show that, for a suitable class of functions a, 
"touch-down" cannot occur in finite time for the transformed problem and as a corollary we obtain 
existence of global solutions to ([l])-(!4l). Furthermore, as a byproduct of our analysis, we also establish 
that the solutions are bounded in uniformly with respect to time. This result applies in particular 
to 

o-i^) = 7-^ ^5 P G (-00, 1] , and a(r) = — — -— , a G (-00, 1] . (5) 

^ ' (l + r)P ^ ^ ^ ' (1 + r)(log(l + r))" ^ ^ ^ ^ 

In the subsequent Section [5] we show that, whatever the value of the mass M is, singularities of 
solutions appear in finite time for some diffusion coefficients between (1 + r)"^ and (1 + r)~^,p > 1, 
namely 

«W<(Y^, Pe(l,2], and a{r)< ^^^^^ (log(l + r)r ' "^^^'^l- (6) 

Thus we eliminate some candidates for the critical nonlinearity (if any!) for ([I])-(jlj). Observe that 
the gap between the functions listed in ([5]) and in ([6]) is quite narrow. The proof of the finite time 
blow-up relies on some virial identities which differ from the ones used originally for the classical 
SP system corresponding to a(r) = 1 P, [I7j. Here again, while the virial identity is written in 
Section in terms on the indefinite integral of u, we mention that the initial guess for this identity 
was made on the transformed equation introduced in Section [21 thus emphasizing once more the 
usefulness of this alternative formulation. Virial identities of the same type have recently allowed 
us to obtain some finite time blow-up results for the parabolic-parabolic quasilinear one-dimensional 
Keller-Segel system [8] (in that case, the equation ([2]) for v is replaced by a parabolic equation). Up 
to our knowledge this is the only result concerning blow-up in finite time for the parabolic-parabolic 
Keller-Segel system besides the remarkable (and more precise) result in [12]. 



2 Change of variables. Blow-up means touch-down. 

In this section we introduce the change of variables that enables us to reformulate the quasilinear 
SP system in such a way that we arrive at a one-dimensional quasilinear parabolic equation 

with an affine right-hand side. It involves the space derivative / of the inverse of the cumulative 
distribution of u. Its surprising applicability to has its origins in the fact that it transforms 

the drift term dr^ {u d^v) in a much simpler one, see f|T5|) below. To the best of our knowledge, such 
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a change of variables has not appeared so far in the hterature in the study of global existence of 
solutions. However, since the g-Wasserstein distance between two functions is related to the L^- 
norm of the difference of the inverses of their indefinite integrals in one space dimension (see, e.g., 
[20] ) , our approach is somehow related to some recent works on degenerate diffusion equations where 
Wasserstein metrics are used to investigate the large time behaviour and the speed of propagation 
of solutions to these equations (see, e.g., [31 HI [TUl [20] and the references therein). 

Let T > and assume that u G C([0,T) x [0, 1]) H C^'2((0,T) x [0, 1]) is a classical solution to 
([T])-(|1D such that n > in [0,T) x (0,1). For t e [0,T), we define the indefinite integral of u{t) 
vanishing at zero by 

U{t,x):= / u{t,z)dz, xe[0, 1], (7) 
Jo 

and observe that the positivity of u warrants that x t— U{t,x) is an increasing function from [0, 1] 
onto [0,M]. For each t G [0,T), let y t— F{t,y) denote its inverse function which is well-defined, 
maps [0, M] onto [0, 1], and satisfies 

U{t, F{t, y)) = y, (t, y) G [0, T) x [0, M] . (8) 

We next set f(t,y) := dyF{t,y) for {t,y) G [0, T) x (0,M). Differentiating ([8]) once with respect to 
y, we first see that 

/(t, y) u{t, F{t, y)) = l, {t, y) G [0, T) x [0, M] . (9) 

Next, differentiating ([9]) with respect to the space variable and differentiating ^ with respect to the 
time variable we arrive at 

d.,u{t, F{t, y)) fit, yf + u{t, F{t, y)) dyf{t, y) = , {t, y) G (0, T) x (0, M) , (10) 

and 

dtU{t, F{t, y)) + u{t, F{t, y)) dtF{t, y) = , (t, y) G (0, T) x (0, M) . (11) 
Now, integrating ([1]) and ([3]) with respect to space over (0, x) we obtain 

dfU {t, x) = a{u{t, F{t, y)) dxu{t, x) — u{t, x) dxv{t, x), —dxv{t, x) = U{t, x) — Mx, 

for {t,x) G (0,T) X (0, 1), so that 

dtU{t, x) = a{u{t, x)) d^u{t, x) + u{t, x) {U{t, x) - Mx) , (t, x) G (0, T) x (0, 1) . (12) 
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Therefore, for {t,y) G (0,T) x (0, M), we have 

dtUit, F{t, y)) = a{u{t, F{t, y)) d^t, F{t, y)) + u{t, F{t, y)) {y - MF{t, y)) (13) 
by dH]). We next use Q, ffTU]l . and ffTTl) along with the positivity of u to deduce that 

(+ ]7(+ w ^ T7(+ \ ^(^) Fj't. y))aiu{t, Fit, y)))dyf{t, y) 
-u{t, F{t, y)) dtF{t, y) = Jifyy ^ ^^^^^ ~ y>' 

and thus 

dtFit, y) = a (y^) dyfit, y) - y + MF{t, y) , (t, y) G (0, T) x (0, M) . 

Differentiating both sides of the previous identity with respect to y and introducing 

^'(r):=^a(-) for any r > , ^(1):=0, (14) 



we reahze that / solves 

dtf = d^^if) -1 + Mf, {t,y)e (0, T) X (0, M) . (15) 
Next, recalling (fTOl) . the boundary conditions (jlj) become 

a,/(t,0)=9,/(t,M) = 0, tG(0,T), (16) 
while the initial value for / is given by (j4]), ([8]), and ([9]), and reads 

/(0,y) = /o(y):= .p^n ?/e(0,M). (17) 

Mo(^(0,|/)) 

Moreover the conservation of mass yields 



fit, y)dy = F{t, M) - F{t, 0) = 1 , t G (0, T) . (18) 

We have thus established the following result: 

Proposition 1. Let T > and consider a classical solution u G C([0, T) x [0, 1]) nC^'^((0, T) x [0, 1]) 
to ([I])-(j4]) such thatu > in [0,T)x(0, 1). Defining U and F by ([7]) and ([8]), respectively, the function 
f G C([0,T) X [0,M]) nC^'2((0,T) X [0,M]) c/eiznec? fty / := dyF is a classical solution to (^-(^ 
satisfying the conservation of mass (ITSl) . 
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Remark 2. Since the main issue of our analysis is whether solutions to (^-(^ exist globally or blow 
up in finite time, we emphasize here that, according to iQ, finite time blow-up atTmax of a solution 
u to corresponds to finite time "touch-down" at T^ax for the related solution f to f|T5l) - f|T7j) . 

namely, 

lim min {fit.y)} = 0. (19) 

In other words, f ceases to he positive at time T^ax o,nd reaches the zero value somewhere in the 
interval [0, M] . 



3 Well-posedness. 

We first recall the following result P, Theorem 1.3]. 

Proposition 3. Consider a positive function a G C^([0, cxd)) and a nonnegative initial condition 
Uq G C([0, 1]) such that (uq) = M for some M > 0. Then there exist a time T^ax ^ (0, C)o] and a 
unique solution 

{u,v)eC{[Q,Tmax) X [O, l] ; R^) H C^'' ( (O, T^ax) X [0,1]; M^), n>o, 

to such that {u{t)) = M for any t G [Q.Tmax)- Moreover, if T^ax < oo then \\u(t)\\^ — > oo 

as t Tmax- By the strong maximum principle we also have u > in {0,Tmax) x (0, 1). 

Proposition [3] does not apply to the case a(r) = 1/r, r > 0, because of the singularity of a at 
zero. In that case, we nevertheless have the following result under the stronger assumption that the 
minimum of Uq is positive. 

Proposition 4. Consider a positive function a G C^((0, oo)) and an initial condition uq G C([0, 1]) 
such that Uq > mo > and (uq) = M for some M > and mo G (0,M). Then there exist a time 
Tmax £ (0, oo] and a unique solution 

{u,v)eCi[0,Tmax) X [0,l];R=^)nCi'2((0,T„,,) x [0,1]; R=^), 

to (HD-dlD such that 

Mm 

= M and u{t, x) > , ° r > for all (t, x) G [0, T^ax) x [0, 1] . (20) 

mo + e^"^{M — mo) 



In addition, T^ax < oo implies that \\u{t)\\^ — > oo as t ^ Tj 



max ■ 
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Proof. We first prove the existence of a solution by a standard approximation argument: more 
precisely, for e G (0, 1) and r > 0, we put a^ir) := a{r + e). As is a positive function in C^([0, oo)), 
we infer from Proposition [3] that there exists a unique classical solution to with 

instead of a which is defined on its maximal existence time interval [0,T^^^). Moreover by standard 
parabolic regularity results [TB] we have 

||-u^||ci+a,2+2«((o,T)x[o,i]) < C{T) ||'U^||L°°((o,r)x(o,i)) (21) 

for some a G (0, 1/2) and T G (0,T^„^.)- By satisfies 

dtu' = (a {u') d^u') - d^v' d.,u' + {u'f - u'M , (t, x) G (0, T^,J x (0, 1) , (22) 
9,.n^(t,0) = 9,n^(t,l) = 0, t G (0,T^,J . (23) 

Introducing 

^^^^ _ J«^U_^ := , i'fr!^ ^ e (0,min{T^„^, ||uo||^^}) , 

1 - ||mo||oo t mo + e*^*(M - mo) 

we notice that S is a supersolution to (!22|) - (!23|) while a is a subsolution to (!22l) - (l23l) with 

a(0) = mo < no(x) < ||no||oo = S(0) , a; G [0, 1] . 
The comparison principle then warrants that 

a{t)<u^{t,x)<m. it,x) G (O,min{r^,,,ho||;o'}) x (0,1). (24) 

According to Proposition [31 a first consequence of fl2^ is that T^^^ > ||^^o||^^- We next infer from 
f ET]) . (^^, and the Arzela-Ascoli theorem that there are a nonnegative function u and a subsequence 
of (m^) (not relabeled) such that 

u'^u in C{[0,\\uo\\^')x[0,l])nC'^mO,\\uo\\^')x[0,l]) (25) 

and 

a(t) <M(t,x) < S(t), (t,x) G (0,||mo||;o') X (0,1). (26) 

If V denotes the unique solution to —d^v = m — M in (0, ||mo||^^) x (0, 1) with homogeneous Neumann 
boundary conditions and mean value zero, it readily follows from (1221) . and that 

(m, w) is a classical solution to with the expected diffusion coefficient a{u). Next, either 
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||'?^(^)||oo — oo as t — >• 1 1 Mo 1 1 and the existence proof is complete, or there is a sequence (ti)i>i 
such that ti \\uo\\^ and (||M(ti) ||oo)i>i is bounded. In the latter, we carry on the construction of 
u starting from u(tio) for a suitable choice of iq > 1. 

Concerning the uniqueness issue, we shall use the Wasserstein distance. Consider an initial 
condition uq fulfilling the assumptions of Proposition H] and let Ui and U2 be two classical solutions 
to enjoying the properties listed in Proposition H] on their respective maximal existence time 

intervals [0,Tmax,i) and [0,Tmax,2)- Thanks to the smoothness and positivity of ui and U2, the 
transformation defined in Section [2] is well-defined for both of them. Consequently, according to 
Proposition [H we may associate to Ui a solution fi to 

dj, = aJ(v|/(/,))-l+M/,, {t,y)e{0,Tma.,^)x{0,M), 
dyMt, 0) = dyf,{t, M) = , t e (0, T^„,,i) , 

/.(o,2/) = My), y&{o,M), 

for i = 1,2, the functions \E' and /o being defined in (fT4l) and Proposition [H respectively. Now, if 
T G (0, min {Tmax,i 5 ^max, 2}) 5 fi and /2 are smooth functions which are bounded from above and 
from below by positive constants (depending on T) and standard arguments ensure that 

||/i(t)-/2(t)||i< ||/i(0)-/2(0)||ie*^* = for t G [0, T] 

(see, e.g., [191 Theorem 11.2]). Consequently, /i = /2 and also Fi = F2 in (0,T) x [0,M], where 
Fi denotes the indefinite integral of fi vanishing at zero as in Section [2J This readily implies that 
ui = U2 in (0,T) X [0, 1] and the proof of Proposition H] is complete. □ 

Remark 5. In the following, we will always assume that there is mo G (0, M) such that Uq > tuq, 
though it is not needed to apply Proposition\^ when a G C^([0, 00)). This is not a restriction since, 
if a & C^{^i 00)), Propositionl^ warrants that u{t) is positive and continuous for each t G (0,Tmax)- 
Instead of Uq, we take as an initial condition uifo) for some t^ G (0,Tmax) which satisfies the above 
positivity requirement. The corresponding solution to is then {t,y) ^ u{to + t,y) for {t,y) G 

{0,Tmax — to) X [0, 1] and Tmax o.nd Tmax — to are simultaneously finite or infinite. 

4 Global existence. 

As mentioned in the introduction we first establish that the diffusion coefficient a{r) = 1/r is not 
critical in the sense that all solutions to are global in that case (Theorem [6]). We in fact prove 
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that the corresponding solution / to the transformed equation (1151) cannot reach zero in finite time: 
to this end, we deduce an L°°-bound on log / from a natural Liapunov functional associated to f|T5|) . 
In a second step, we extend our analysis to a wider class of diffusion coefficients a (Theorem [7]). 

Theorem 6. Assume that a{r) = 1/r and consider a non-negative function uq G C([0, 1]) such that 
uq > rnQ > and (uq) = M for some M > and ruQ G (0,M). If{u,v) denotes the corresponding 
solution to then T^ax = oo, the trajectory {{u{t),v{t)) : t > 0} is bounded in L°°{0, 

and there is fi > such that u{t, x) > fi for all {t, x) G [0, oo) x [0, 1]. 

Global existence is actually true for a wider class of nonlinear diffusion coefficients a such as 

O'ir) = -, ^ , V G (— oo, 11 , and a(r) = — ; -— , a G (— cxd, 11 . 

(l + r)p' ^ ' ^' ^' (1 +r)(log(l +r))" ' ^ ' ^ 

Indeed, the above examples fulfil the assumptions of the following result: 

Theorem 7. Assume that a G C^((0, oo)) is a positive function such that a ^ L^{1, oo) and, for each 
e G (0, oo), there is > for which 

a{r) < e ra{r) H — - for r G (0, 1) . (27) 
r 

Consider an initial condition uq G C([0, 1]) such that uq > mo > and (mq) = M for some M > 
and tuq G (0, M). If {u,v) denotes the corresponding solution to ([I])-(j4]) given by Proposition 
then Tfnax = oo and the trajectory {{u{t),v(t)) : t > 0} is bounded in L°°(0, 1;]R^). In addition, if 
a ^ -^^^(0, 1), there is fi > such that u{t, x) > fi for all (t, x) G [0, oo) x [0, 1]. 

Though the proofs of Theorems [6] and [7| follows the same steps, the generality of the latter brings 
additional technicalities in the proof. We will thus first present the proof of Theorem [6] to illustrate 
the main ideas and then proceed with the proof of Theorem [71 

Proof of Theorem [6l According to Proposition HI it suffices to provide a control on the L°°-norm 
of u on any finite time interval to ensure global existence. By Proposition [1] and Remark [2] such a 
control can be obtained by showing that the corresponding solution / to (fT5!) - (fT7j) is bounded away 
from zero on any finite time interval. To achieve the latter, we need some more steps. Recalling that 
\E'(r) = logr in that case, we start with the following observation. 

Lemma 8. The function 

^^^y-= I (^l\dy\ogf{t,y)\^ + \ogf{t,y)^ dy, t e [0,Tmax) , (28) 
is a non-increasing function of time. 
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Proof of Lemma [5]. Multiplying (fT^ by 9^ log / and integrating with respect to space over (0, M) 
we arrive at 

flty) "^^^"/o ^^''Sf(t^y)9tdylogf{t,y)dy-J^ d^log f{t,y) dy+M dj{t,y)dy. 

Owing to the conservation of mass ( |T8i) . the last term of the right-hand side of the above identity 
vanishes and we end up with 

dr 

^W + l f{t,y)mogf{t,y)\'dy = 0, (29) 

from which Lemma M readily follows. □ 

At this point, we notice that, if a bound from below turns out to be available for L, we might 
expect to deduce from it an estimate on log/ in L°°{0,T] H^{0, M)). Thanks to the continuous 
embedding of H^{0, M) in L°°(0, M), such an estimate will in turn provide an L°°-estimate on log/ 
from which clearly follows that / is bounded away from zero. The next step in the proof of Theorem [6] 
thus requires a detailed study of L. To this end, we define the functional 

E{g) := J^'' (^^\dygiy)\' + giyij dy for geH\0,M), (30) 

so that Lemma [8] states that t i — > E[log f{t)) is a non-increasing function of time. We next proceed 
to find a lower bound on E. 

Proposition 9. If g e H\0,M) satisfies (e^) = 1/M then 

E{g) >^\\dyg\\l- M log M - , (31) 

and 

\\g\\i < 2 + MlogM + M^/^\\dyg\\2 . (32) 

Proof of Proposition [HI We consider g E H^{0,M) such that (e^) = 1/M. Thanks to the 
continuous embedding of H^{0,M) in C([0,M]), belongs to C([0,M]) and the constraint (e^) = 
1/M and the mean value theorem ensure that there exists Yg G [0, M] such that e^'-^"-' = 1/M. Then 

I'M I'M 

/ 9iy)dy = / (giy) - giYg)) dy + Mg{Yg) 
Jo Jo 

I'M py 

> dyg{z) dzdy - M log M 

Jo JYg 

> -M^/^\\dygh- MlogM. (33) 
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Thanks to fl33p we can estimate E{g) from below and obtain 

Eig) > ^\\dyg\\l + ^\\dyg\\l-M'/'\\dyg\\2 + M'-M'-M\ogM 
> \ \\dy9\\l + \ {\\dygh-2M'/'Y -M'-M log M, 
hence (IMl) . It also follows from ([3SD that 

\\g\\i = 2/ ma.x{0, g{y)} dy - g{y) dy 
Jo Jo 

pM 

< 2 / (^e^^-{o,9{y)} dy + M log M + M^I'^Wdygh 

Jo 

I'M 

< 2 e'^ (y) dy + M log M + M'/^dygh, 

Jo 

from which ([32]) follows as (e^) = 1/M. □ 
We now complete the proof of Theorem [6] by noticing that Proposition [9] implies that log / G 
L°°{0,T] H^{0, M)) for any T G {0,Tmax) with a bound that depends neither on T nor on T^ax- 
Indeed, since (e'°s^) = 1/M by ([I8]), we infer from ([28]) and ([3l]) that, for t G [0,T„a^), 

^||a,log/(t)||^ < Eilog fit)) + M log M + M'< Eilogfo) + MlogM + Ml 

In view of f[32]) . the previous bound entails that log / belongs to L°°{0, Tmax', -^^(0, M)) and we finally 
conclude with the help of the Poincare-Wirtinger inequality that 

||log/(t)||^,i(o,M)<C(M)(l + E(log/o)) for te[0,T^ax). (34) 

Owing to the continuous embedding of H^{0, M) in L°°(0, M) and ( 134]) . we conclude that there are 
constants < a < /3 such that a < fit,y) < [3 for {t,y) G [0,Tmaa) x [0,M]. This property clearly 
warrants that Tmax = oo and ends the proof of Theorem [6] for a(r) = 1/r. □ 

Proof of Theorem [Tj. We follow the same strategy as in the previous case and aim at showing 
that a similar argument prevents the solution / to (fT5]) - (fT7)) to hit zero in finite time, the function 
being given by ( fT4l) (see Section [2]). Owing to the positivity and the non-integrability over (1, oo) 
of a, we note that \E' is an increasing function from (0, oo) onto its range with inverse \E'~^ and 

\1/ maps (0,1) onto (— oo,0) with lim \l/(r) = — oo . (35) 

We next report the analogue of Lemma [HI 
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Lemma 10. The function 

-M I'M 



-1 I'M rlVl 

U{t):=-j^ \dy^ifit,y))f dy + {m{f{t,y))-M^,{f{t,y)))dy 

is a non-increasing function of time on [0,Tm,ax), the function being defined by 

^i(l):=0 and (r) := r^'(r) = - a ( - ) , r G (0, cx)) . (36) 

r \r J 

Proof of Lemma IIOI We multiply f|T5l) by dt'^lf) and integrate with respect to space over (0, M) 
to obtain 

I'M pM 

/ vl>'(/) dy = - dy<il{f)dtdy^{f) dy 

Jo Jo 

I'M I'M 

- dt<^U)dy + M fm\f)dtfdy, 
Jo Jo 



hence 

dU 



{t)+ / ^'{f) \d,f\'dy = Q. 
Jo 



dt 

and Lemma fTOl follows in view of the positivity of □ 
Observe next that, owing to ( !27l) with e = 1/M, we have 



M ^ s = — a - < — 
s \s J s 



(ilk " (s) ^ ^^^^^ ^) - ^ ^ ^ ^ 



so that, after integration over (1,t), r > 1, M\l/i(r) < \l'(r) + MnijM r. Consequently, since < 
on (0, 1], we have 

I'M I'M 

-M / m,{f){t,y)dy > -M / y)) v[/i(/(t, y)) rfy 

Jo Jo 

I'M 

> - / l[i,oo)(/(t,y)) (^(/(t,y)) + MKi/M dy 







"M 

> - I l[i,oo){f{t,y))^ifit,y))dy-MK,/M 







where we have used f|T8l) to obtain the last inequality. Now, recalling fl35|) . we find 



M I'M 
|2 



Li{t) > I- I \dy^{f{t,y))\' dy+ I l^o,i){f{t,y))^{f{t,y))dy-MKyM 
Jo Jo 

^ I'M rM 



> \dy^{f{t,y))\' dy + j l^_^^,0{f{t,y)))^>{f{t,y))dy-MKyM 

> E,{^{fm-M ^i/M, (37) 
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with 



1 

Ei{h):=-\\dyh\\l + l(^^^^){h{y)) h{y) dy , heH\0,M). 
The next step is then to find a lower bound on Ei. 
Lemma 11. For M > 0, we have 



Ei{h) > - \\dyh\\l-M^-M 



* — 

M 



and 



for every h G -f^"'^(0, M) satisfying 



i<M'/'\\dyhh + M 



^ — 

M 



M 



^~\hKy) dy=l. 



(38) 
(39) 
(40) 



Proof of Lemma llll Consider h G H^{0,M) enjoying the property fHU]) . Arguing as in the proof 
of Proposition [9] , we infer from (140!) and the mean value theorem that there is at least one point 
Yh G [0,M] (depending on h) such that ^'~^(/i(Y'h)) = 1/M. Consequently, 



M 



min {0, h{y)} dy 



M 



> 



(min {0, h{y)} - min {0, h(Yh)}) dy + Mmm {0, h{Yh)} 

AI ny 



'0 JYn 

> -M^l^\\dyh\\2~ M 



{-oo,o){h{z)) dyh{z) dzdy - M 



1 

M 



1 

M 



(41) 



Thanks to (14T]) we can estimate Ei{h) from below and obtain 

E,{h) > 1 \\dyh\\l + \ {Wdyhh - 2M'/'Y -M'-M 
hence (l38l) . We next proceed as in the proof of (14T!) and compute 



1 

M 



M 



\h{y)\ dy 



M 



< 



{\h{y)\ - \h{Y,)\) dy + M\h{Y,)\ 

AI i-y 

sigii{h{z)) dyh{z) dzdy + M 



M 



10 JYh 

< M^'^\\dyh\\2 + M 



1 

M 
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thus establishing (15^ . 



□ 



Owing to f|T8l) . fulfils the property fHOl) and we infer from Lemma [TOl fl37j) . and fl38|) that, 
for t G [0,T^ax), 



Thanks to the Poincare-Wirtinger inequality the above two estimates imply that : t G 

[0, Tmax)} is bounded in L°°(0, Tm^x; -f^^(0, M)). Owing to the continuous embedding of H^{0, M) in 
C([0,M]), we conclude that there is K > such that -K < < K for (t, y) G [0,Tmax) x 

[0, M]. Since \& is an increasing function, we actually have \['~^(— Ji) < f{t,y) < \E'~^(A') for 
{t,y) G [0,Tmax) X [0, M] and ^-^(-if) > by ([27]). Recalling Remark [21 we thus have shown that 
Tmax = OO. Finally, since the lower bound on / do not depend on time, the boundedness of the 
trajectory {{u(t),v{t)) : t > 0} in L°°(0, 1;M^) readily follows. Under the additional assumption 
a ^ L^{0, 1), the function \1/ maps (1, oo) onto (0, oo) and the time-independent upper bound on / 
provides a positive lower bound on u thanks to Q. □ 

5 Finite time blow-up. 

In this section, we return to the question of finite time blow-up and show that it can be proved for a 
large class of diffusion coefficients a with the help of a suitable virial identity. Besides providing finite 
time blow-up results for functions a which were not included in [5], it also provides an alternative 
and simpler proof of [H Theorem 3.2]. 

Theorem 12. Let a G C^((0, oo)) be a positive function such that a G L-'^(l,oo) and there is a 
concave function B for which 



\\dy^{fm\l < 4Ei(vl/(/(t)))+4M3 + 4M VI/ 
< 4Li(t)+4MKi/M + C(a,M) 




< 4Li(0) + C(a,M). 
Moreover, combining the previous inequality with fl39|) . we obtain that 

||^(/(t))||i<M3/2 {4L^{0) + C{a,M)f' + M f-J-") 




t G [0, TfYiax^ ■ 




r G (0, oo) 



(42) 



lim 



0. 



(43) 



■oo 



r 
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Then, for any M > 0, there is at least an initial condition Uq G C([0, 1]) satisfying Uq > tuq > 
and (uq) = M for some mo € (0, M) for which the first component u of the corresponding solution 
to ([I])-(jl]) blows up in finite time, i.e. T^ax < oo. 

Observe that Theorem [T2l applies when a(r) < C(l + r)"^ for some C > and p G (1,2] (with 
B{r) = Cr{l + rY~^ /{p — 1)). In that case, we thus recover the blow-up results from [9] by a 
different method. However, Theorem [12] also applies if there exist C > and a G (1,2] such that 
a{r) < C(l + r)-i (log(l + r))~" (with B{r) = Cr(log (1 + r))i~"/(a - 1)). 

Proof of Theorem [lH Let g > 2. Introducing Lg{t) := \\U{t)\\yq (recall that U{t) is the 
indefinite integral of u{t) vanishing at zero, see Section [2]), we infer from fll2p . (1421) . the fact that 
U{t, 0) = M - U{t, 1) = 0, and the negativity of A that 

dLn 



dt 



[ U^-^ d^A{u) dx+ [ u{U - Mx) U"'^ dx 
Jo Jo 



[U''-' A{u)]l~l - {q - 1) I U"~^uA{u)dx + 



q + l 



x=l 



x=0 



M 



X 



x=l 



+ M L, 



J x=0 



< {q 



U'^-^ Biu) dx + M L„ 



By Jensen's inequality with measure B{u) dx, we have 



(44) 



U''-^ Biu) dx 







Biu) dz 







< 



1 \ 2/g 

Biu) dx 



'(f;9)('/-2)/g B{u) dx 
B{u) dz 

1 \ {9-2)/g 

U'' Biu) dx 



< {qL, 



\(9-2)/<i, 



1 \ 2/g 

B{u) dx j 



\iu) "^"^ 







qL^ 



We now use the concavity of B, the Jensen inequality (with measure dx and If^ dx), and the 
conservation of mass (l20l) to obtain 



n2/g r 



< (gi:,)^'^"')/'' [BiM)]'/" 



B 



B 



j'^W-^B{u)dx < (gL,)^^-')/'' B(^J\dx^ " ^ ^ ^) 



('7-2)/<? 



qiq + i)L, 



qLq, 

(g-2)/g 
(9-2)/g 



('/-2)/q 



q{q + ^)L, 
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where (3{r) = B{r)/r for r > 0. Recalling fl44p and introducing 



AM(r) :=Mr + (g-l) [B{M)] 



for r > 0, we conclude that 



2/q 



(g-2)/g r 



g(g + l)r 



{9-2)/g 



~dt 



(t) < AM{Lg{t)) , te[0 



(45) 



Assume now for contradiction that Tmax = oo. Since Am(0) < by (H3l) . there is i^Af > such that 
AAf(r) < for r G [O,-;?!/]- Then, if Lg{0) < -i^^/, we readily infer from fj45l) and the negativity of A^/ 
in [Oj^^m] that ivg(t) < '&m and dLq{t)/dt < sup^^^j,^] {Aj\/} < for all t > 0. This property implies 
that Lq becomes negative at a finite time and thus contradicts the non-negativity of Lq. Therefore, 
Tmax < oo provided that Lq{0) < ^m- The latter condition is fulfilled as soon as Uq is sufficiently 
concentrated near x = 1, whatever the value of the mass M is, and the proof of Theorem [12] is 
complete. □ 
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